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ABSTRACT 

We present a mechanism for the formation of the low frequency 1 // magnetic spectrum based on numerical 
solutions of a shell reduced-MHD model of the turbulent dynamics inside the sub-Alfvenic solar wind. We 
assign reasonably realistic profiles to the wind speed and the density along the radial direction, and a radial 
magnetic field. Alfven waves of short periodicity (600 s) are injected at the base of the chromosphere, pen- 
etrate into the corona and are partially reflected, thus triggering a turbulent cascade. The cascade is strong 
for the reflected wave while it is weak for the outward propagating waves. Reflection at the transition region 
recycles the strong turbulent spectrum into the outward weak spectrum, which is advected beyond the Alfvenic 
critical point without substantial evolution. There, the magnetic field has a perpendicular power-law spectrum 
with slope close to the Kolmogorov -5/3. The parallel spectrum is inherited from the frequency spectrum of 
large (perpendicular) eddies. The shape is a double power-law with slopes of ^ -1 and -2 at low and high 
frequencies respectively, the position of the break depending on the injected spectrum. We suggest that the 
double power-law spectrum measured by Helios at 0.3 AU, where the average magnetic field is not aligned 
with the radial (contrary to our assumptions) results from the combination of such different spectral slopes. At 
low frequency the parallel spectrum dominates with its characteristic 1 // shape, while at higher frequencies its 
steep spectral slope (-2) is masked by the more energetic perpendicular spectrum (slope -5/3). 
Subject headings: Magnetohydrodynamics (MHD) — plasmas — turbulence — solar wind 



1. INTRODUCTION 

At heliocentric distances of about 0.3 AU the magnetic 
field spectrum in the fast streams of the solar wind has a 
form of a double power-law with a break at about / - 
5 10"^ Hz and a slope of approximatel y -1 and -5/3 at lower 
and higher frequencies respectively (iBava ssano et al.l 119821 : 
iDenskat & Neubaueii[T98l iBruno & CarbonelllOOsT Fluc- 
tuations in these two ranges have a different evolution with 
distance (Bavassano et al. 1982; Marsch & Tu 1990; Roberts 
[1992). The low frequency part approximately follows the 
WKB behavior, E oc r"^, dictated by the solar wind expan- 
sion. The high frequency part instead decreases much sharply, 
its energy content being depleted by the turbulent cascade, 
and maintains approximately the same spectral slope. There- 
fore, the frequency break shifts to lower and lower frequency 
as the heliocentric distance increases, the 1// spectrum oc- 
cupying the range 3 K)"^ Hz < / < 8 lO^^ Hz at 1 AU 
(iMatthaeus & Goldsteinlll986h . The energy in the low fre- 
quency tail progressively contributes to the turbulent heating 
of the solar wind. The spectral evolution can be understood 
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as a competition between the expansion timescale and the cas- 
cade timescale, regulatin g the decay of e nergy respectively at 
small and at large scales ( Tu et al .11 19 841) . but the origin of the 
1 // spectrum is still a matter of debate. 

There are some indications that it may have a genuine solar 
origin, reflecting the emergence, cancellation, and sinking of 
the magnetic field at the photosphere: spectra built from mag- 
netogram data at low-intermediate latitudes show a 1 // slope 
at low wavenumbers, which also correlate to the spectrum 
of de nsity fluctuations in the heliosphere (^ Matthaeus et al.l 
l2007h . In this view, the formation mechanism relies on mag- 
netic reconnection, occurring over a hierarchy of scales as a 
stochastic process with some self-similar properties. Recon- 
nection, this time in the corona, has also been invoked as the 
underlying process, since the timescale associated to the re- 
structuring of coronal magnetic field is about 1 day , i.e. in the 
range of observed frequencies (IClose et al.l 12004b . Recently 
the formation of a 1 // spectrum has been observed in homo- 
geneous incompressible MHD simulations of very long du- 
ration, probably originating from an i nverse cascade not as- 
socia t ed with well defined invariants (iDmitruk & MatthaeusI 
120071; IDmitruk et al.l 120 lib . Other ideas rely on large scale 
properties of the corona and solar wind: reflection was sug- 
gested as a possible mechanism to obtain a fc"' spectrum from 



an isotropic cascade of Alfven wa ves waves comin g from the 
Sun in the expanding solar wind dVelliet al.lll989l) . The lin- 
ear evolution of Alfvenic pulses excited in the corona leads 
to a signal (outside the Alfvenic critical point) which has 
periodicities of about 15-30 minutes (ringing of the corona, 
iHollweg & Isenberg 2007). Though this period appears too 
low to account for the entire range of observed frequencies, 
as suggested by the authors, the ringing of the corona could 
play a role in the formation of the 1 // spectrum. 

In this Letter, we combine the latter two ideas, the ring- 
ing of the corona and nonlinear interactions with reflected 
waves, to study the formation of the magnetic field spectrum 
advected by the sol ar wind. To this aim , we use a Shell Re- 
duced MHD model dVerdini et alj2009l) to follow the onset of 
turbulence resulting from the coupling between Alfven waves 
propagating in the chromosphere, corona, and sub-Alfvenic 
solar wind. The solar wind profiles (Alfven speed, velocity 
and density profiles) are given steady nurn erical solutions of 
the slow wind model by iPinto et al.l (120091) . Such slow wind 
solutions are chosen to illustrate the mechanism responsible 
for the formation of the 1 // spectrum, fast wind solutions will 
be considered in a forthcoming paper 

In homogeneous turbulence with equal energy in both 
Alfven species the forcing controls the development of a weak 
(or strong) cascade, depending on whether the parallel Alfven 
time is smaller (or not) than the nonlinear time. The energy 
dens ity scales with pe rpendicular wave number as k]^ or fc^ 
(e.g. IVerdini & Grap pin 2012) for both species. However, in 
the inhomogeneous stratified open corona and solar wind, the 
two Alfven wave species can in principle have very different 
amplitudes. We will see that they follow different cascade 
regimes (weak/strong) at the same time. This is the key of 
the results obtained in the present Letter We show that for 
realistic conditions a 1 // magnetic spectrum forms as a con- 
sequence of the strong turbulent cascade of reflected/inward 
propagating waves, which reflect back (ringing) at the transi- 
tion region (TR) and are advected outside the Alfvenic criti- 
cal point without substantial modification, since the outward 
propagating waves experience a weak cascade. 



2. EQUATIONS AND PARAMETERS 

The model equations are obtained from the MHD equations 
by separating the l arge-scale stationary fields and small-scal e 
fluctuating fields (iHeinemann & Olberti 119801: IVeflil Il993h . 
Among the large-scale fields, the wind speed (t/) and density 

i p) are the s olution of a ID Hydro -dynamic solar wind model 
Pinto et all 120091: iGrappin et al] 12010.) with specified heat- 
ing functions and assigned flux tube expansion (A - Tiry^). 
Combining p with the magnetic field strength B© at the base 
of the fluxtube, and requiring BA - const one obtains the 
Alfven speed {Va - B/ -sjA-np). The profiles of U, Va, p, and 
;F depend only on the radial distance r - R/Rq and are plotted 
in the left panels of Figure [T] the bottom and top boundaries 
are at r^c = 1.004 and r,op = 19, the Alfvenic critical point 
is found at r^ ~ 17 by choosing B© = 10 G. The small-scale 
incompressible velocity (u) and magnetic (b) fluctuations are 
orthogonal to the radial B and can be expressed through the 
usual Elsasser variables z^ - u + b/ -\j4-np. We further sim- 
plify the equations by replacing the nonlinear terms and the 
pressure term at each position by a 2D shell model, so that the 



final equations read: 

dt dr 4 \p drj 



1 1 1 dp I dA\ ^ , -, , 



where we take equal viscosity and resistivity (v) and r,f 



de- 



notes the nonlinear interactions. The index n labels modes 
having perpendicular wavenumber k„ - ko2'\ which de- 
fine the radius of concentric shells filling the (perpendicular) 
Fourier space. The largest transverse scale follows the flux 
tube expansion ^o('") = ^oo/ VS(r), wavevectors are given 
by kj_ = k„(r) = k{)(r)2", so that our Fourier space shrinks 
with increasing r. Two complex fields are assigned to each 
shell, z^(r, t) = z^(k„, r, t). They have the dimension of a ve- 
locity and kjJ'l^/2 are the respective energies per unit mass in 
the shell n. We recall that nonlinear interactions are local in 
Fourier space, T„pq oc 'Lp^knZpZa, with p ~ q ~ n (the ex- 
plicit expre ssion and coefficients for t he 2D shell model may 
be found in lGiuliani & Carbonelll998h . 

Open boundaries are used at r^o, and r,op. The free param- 
eters are the input wave amplitude Zg, the perpendicular in- 
jection scales koQ (the forcing being imposed on 1,2,4 A:o0), 
and the correlation time of the forcing signal Tf. We choose 
a strong forcing, i.e., Tf ^ f^^ by assigning: zj = 10 km/s, 
koQ = 2;7r/ 34000 km"', and T/ = 600 s. With these parameters 

C = i/^%-500s-ry. 

Decreasing the free parameter Bq shifts the profile of Va in 
Figure [T] to lower values, thus decreasing r^ and resizing the 
regions where one linear term dominates over the other linear 
terms (propagation, WKB, and reflection, respectively the II, 
III, and IV terms on the LHS of Eq. IB- 
Starting from the solutions zfj{r,t) of Eq. [T]we define fre- 
quency and wave number spectra. We denote by z^{r, f) their 
Fourier transform with respect to time. From this we define 
successively the spectra Ef^if) - E'^in,/) and the associated 
reduced spectra at each radial distance r. 

EtiD^Kirjtlk^ (2) 

E}(f)^jE^dk^ (3) 

E^JkJ^JE^df (4) 

where we have omitted the dependence on r. Equation ^ is 
the frequency spectrum of a given perpendicular mode which 
yields the total energy E'^(r) - 2„ \z'l;\^ - J E^dk±df. Equa- 
tions (l3]|4| define reduced spectra, which accordingly yield 
J E'^df - J E'^dk^ - E^ (hereafter we drop the subscripts 
/, k^ when the dependence is explicit). 

The correlation time t^^^iir, k^) is defined as the full- 
width-half-maximum (FWHM) of the autocorrelation func- 
tion ACiz"^) computed at each position r, while the nonlinear 
time f^^(r, kj_) is the time av erage of the eddy turnove r times 
built from the Elsasser fields (iDobrowolny et al.1ll980b : 



t-„,(kJ=FWUM[AC(zn,] 



''cor\ 

t%LikA.)^{ll{kAz^ik^)\)) 



(5) 
(6) 



Again we omitted the dependence on r, ACiz) - J z(t)z(t - 
t')dt' is the autocorrelation function and (...) stands for a 
time average. The amplitudes of the fluctuations are com- 
puted by integrating along the perpendicular wavenumbers 
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Fig. 1 . — Left panels: solar wind model. Profile.s of the Alfven and wind 
speed (top) and of the density and over-expansion 'F = A(r)/r- (bottom). 
Right panels: snapshots of the amplitude of the Elsasser variables (top) and 
of the velocity and magnetic fluctuations (bottom). The Alfven and wind 
speed are also overplotted as dotted lines. 
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Fig. 2. — Top. Perpendicular spectra ^^(^j^) at r,gp. The input spectrum 
E* at r],Q, is also shown (dashed line). Bottom. Frequency spectra £* at r,„p 
for perpendicular wavenumbers <:o2" (with n increasing from top to bottom). 
Symbols mark the CB width A^,g (see text). See the electronic edition of the 
Journal for a color version of this figure. 



z^ - -\/-^„|z„(r, f)P- The same definitions hold for u and b, 
the latter will be expressed in velocity units from now on (i.e. 
b ^ b/ sjA-np). The equality of the two timescales defines the 
so called critical balance (CB) condition: 



fcorW ~ ^NlW^ 



(7) 



that is supposed to hold for strong turbulence. Since the 
correlation time is the inverse of the width of the frequency 
spectrum, we can define an equivalent CB width A;!;^ = 1 /f^^. 
According to anisotropic turbulence theories 
(iGoldreich & Sridha^ 119951) . this width constrains the 
frequency spectrum, since at a given perpendicular scale kj_ 
there is little energy for / > Acb, while spectra are flat for 
/ < Acb. 

3. RESULTS 

Figure[T]shows a snapshot of the amplitude of Alfven waves 
(top right) and of the kinetic and magnetic fluctuations (bot- 
tom right). In the corona z*' - 1. 5b, b - 1.5m. The pro- 
file of z^ follows the wind speed profile (background dotted 
line) and is much larger the inward/reflected wave z". z^ is 
a smooth function of distance, while z^ displays small scale 
parallel structure, seen also in the radial profiles of u and b. 

The spectral densities vs perpendicular wave number and 
frequency at r,op = 19 are summarized in Figure |2] In the 
top panel the time-averaged perpendicular spectra E^^ik^} are 
plotted along with the input spectrum at the base of the chro- 
mosphere (£(|). Both perpendicular spectra show a well de- 
veloped power-law that extends about two decades with a 
slope -5/3. The ratio E^/E]_ is about constant in the iner- 
tial range. The input perpendicular spectrum E^ (dashed line) 
appears in the middle of the inertial range of the spectra at 
r,op, as the flux tube expansion has strongly expanded the per- 
pendicular scales between the surface and r,op (by a factor 
100). 

In the bottom panels of Figure |2] the frequency spectra, 
E^if), are plotted for each perpendicular mode n. Symbols 
mark the equivalent CB width A^^. The E' spectra (right) are 
very well bounded by the CB condition, falling off' sharply in 
the weak turbulence range (/ > A^^), and having flat spectra 



in the strong turbulence range (/ < A^^^), with few exceptions 
for large-scale eddies (n = 0, 1,2). On the contrary E'^ has 
a substantial over-excitation of high frequencies at large per- 
pendicular scales and in the whole weak range, which reflects 
the fine-scale (parallel) spatial structure of z^ seen in Figure[T] 
This over-excitation has an approximate slope 1// and ex- 
tends from the CB boundary to about l/Tf ^ 2 lO"-' Hz. 

The slope of the frequency spectrum has been found to 
depend on turbulence forcing strength Tfl^l,, in homoge- 
nous shell RMHD simulations ( Verdini & Grappin 2012), and 
to extend beyond forced scales. The situation is different 
here because we cannot really control the forcing turbulence 
strength: the nonlinear time for z^ depends on the reflected 
amplitude z", an output of the simulation. An a posteriori 
estimate yields a turbulence strength t'^grl^NL < 1 so that z^ 
is indirectly forced in the weak regime. Coming back to the 
strong turbulence range, one can see that large eddies have 
an approximate slope of 1// in both E^ and E^: these low- 
frequency fluctuations are lineariy coupled by density gradi- 
ents (reflection) that force the spectra to resemble each other 

The turbulence regime at different heliocentric distances 
can be identified in Figure [3] where we plot the corre- 
lation time vs the nonlinear time for three wavenumbers 
(2.k(), 8^0, 128feo) at heliocentric distances, running from left 
(r = 2) to right (r = 19) along the solid lines. The critical 
balance condition Eq. |7] is drawn as a dotted line and sepa- 
rates the strong and weak regimes (above and below respec- 
tively). The reflected wave z" is always in a strong turbulence 
regime and follows the CB condition at any position, except 
for a small offset at low wavenumbers due to the linear cou- 
pling (reflection) at large scale. On the contrary the correla- 
tion time of z^ is almost horizontal, i.e. almost independent 
of the nonlinear time and of distance, and lies entirely in the 
weak turbulence regime. Since tcr is the inverse of the width 
of the frequency spectrum, one can conclude that the E^ spec- 
trum does not change much as it propagates outward, only 
becoming a bit narrower at large scales. Instead, the E^ spec- 
trum widens while propagating backward from r^, according 
to CB. Once it arrives at the TR it experiences strong reflec- 
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Fig. 3. — Correlation time /*^^ versus nonlinear time f*^ for different helio- 
centric distance (running from 2i?o to 19Rq from left to right on each curve) 
at three perpendicular wavenumbers 2k^, Sir*}^, llHk^j^. The dotted line is the 
CB condition Eq.[7] 
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Fig. 4. — Reduced frequency spectra for £"• at r,op (solid gray and black 
line respectively) and input spectrum E^ at rtoi (dashed line). See the elec- 
tronic edition of the Journal for a color version of this figure. 



tion and feeds the E'^ spectrum, this can be argued by noting 
that t^g,. < Q^ at r = 2 (the low end of the curves). At larger 
distances f^t^^ << f^T^,., showing that E^ is wider than £■„ in 
the whole corona, the difference owing to its high-frequency 
over-excitation. 

The 1 // spectrum appearing at large perpendicular scales 
in E'^ is thus made up of two parts, each occupying about one 
decade from / ^ 10"^ Hz, which originate from two different 
mechanisms: linear coupling with E' at low frequencies (the 
strong turbulence range), and weak turbulent cascade at inter- 
mediate frequencies. This is summarized in Figure S] where 
the reduced frequency spectra £"''(/) of kinetic and magnetic 
energy at r = r,op are shown, (along with the input spectrum 
at the base of the chromosphere E^). The spectrum of £* 
shows the 1 // slope inherited from the large scale eddies in 
£"^, in which one can recognize the above two components, 
and a break appears at about / ^ 10"'' Hz where the slope 
switches to -2. The break coincides with the width of the 
forcing spectrum (1/7/^2 10"^ Hz) and the -2 slope is con- 
sistent with the fall-off at high frequencies in £* ^ E^. Note 
that linear propagation alone would lead to a magnetic spec- 
trum with slope /"'''^ (not shown), thus nonlinear interactions 
are fundamental to obtain the 1// slope. The spectrum of E" 
is practically identical to £* for / > lO"'* Hz due to the same 
over-excitation in the weak regime. For lower frequencies in- 
stead the slope of E" is flatter, having a value of about -1/2. 

4. DISCUSSION 

Why the 1// extends down to low frequencies only in 
E'' can be understood by examining the wave reflection at 
the TR, where density gradients are higher For frequencies 
/ < max[\dVa/dr\] ^ 2 Hz, one can write z" - -(1 - e)z^ 
where e = v^'"''""/V^°'' « 1 is the Alf ven speed con- 
trast between the chromosphere and corona (Hol lwegllI984l 
Ivan Bal legooijen et al. 2011; Verdini et al. 2012). At die TR, 
b ^ z"^ and u « z^, thus reflection transfers the low- 
frequency large-eddy spectrum of E^ to £* but not to E". In 
reality reflection occurs in the whole low corona (for r < 3 Rq 
at frequencies / < 5 10"^ Hz) and is a continuous pro- 
cess. Moreover, reflected waves are made of two compo- 



nents: a classical one, propagating backward in the region 
U < Va and subject to ringing in the low corona, and an 
anomalous comp onent that instead prop agates outward with 
the mother wave dyeUi et al.lll989:.Hollweg & Isenberg|2007t 
I Verdini et al.ll2009h . Nonlinear interactions are not Umited to 
the duration of the encounter of colliding z=^ wavepackets, but 
part of them occur on the common outward propagation path. 
This is why linear and also nonlinear couplings affect the low 
tail of the 1 // spectrum. 

Bearing in mind that z" is entirely generated by reflection 
inside the solar atmosphere, we examine now what happens if 
we increase the input turbulence strength (i.e the control pa- 
rameter Tflt^^j) by either increasing the correlation time Tj 

or decreasing the input nonlinear time f^^ - l/kaQZ^- By im- 
posing a shorter input nonlinear time, no matter if through Zg 
or ^oo, we increase the strength of the cascade for z". The 
strength of the z^ cascade will increase only slightly, since 
the reflected wave will be damped more strongly, yielding 
approximately the same t^^. The spectrum of E"^ will thus 
be affected only slightly, except that now nonlinear coupling 
will dominate over reflection in the low corona, eroding the 
very low frequency tail of the 1 // spectrum in E^. If instead 
we choose a longer correlation time at input, we increase the 
amount of reflected waves thus also increasing the strength 
of the cascade for z^- This time the high-frequency tail of the 
1 // spectrum will be eroded, since a stronger cascade reduces 
the energy in the over-excited weak-turbulence regime of E'^ . 
Finally, one can vary the importance of linear/nonlinear cou- 
pling through the other free parameter Bq. Taking a smaller 
chromospheric magnetic field will lower the Alfvenic critical 
point, reducing the importance of linear terms compared to 
the nonlinear ones in Eq. [1] The changes in the E* spectrum 
are similar to those discussed above for a shorter t^^. Note 
however that halving Bq moves r^ to a distance of 8 Rq, so 
that one is limited to small variations of this parameter A last 
remark concerns the non-local couplings in the Fourier space 
that are neglected in the shell model employed in this work. 
In principle they could change the perpendicular spectra and 
the strength of the cascade. However, our understanding of 
the process of formation of the 1 // spectrum is that it comes 
from the different nature of the two direct z" and z^ cascades. 



not because of an inverse cascade as in the dynamo process. 
We expect these direct cascades not to change when including 
nonlocal interactions. 

The formation of the 1 // spectrum by nonlinear and linear 
coupling in the sub-Alfvenic solar wind appears to be quite 
solid, however its relation to the interplanetary spectrum 
observed at 0.3 AU is not straightforward. Beyond rA the 
solar wind expansion causes the rotation of the magnetic 
field, a slower production of reflected waves, and a decrease 
of wave-numbers perpendicular to the radial direction 
(jGrappin et al. 1993). The latter two actually suggest a 
freezing of the spectra. On the contrary the magnetic field 
rotation causes instabilities and wave coupling/decay that 
could modify the spectrum. Numerical simulations of 
MHD equations in ID and 2D incorporating the effect of an 
expan ding solar wind (Grappin et al. 1993; Grappin & Velli 
11996) suggest that the decay and instability rates are delayed 
with respect to the non-expanding case, so that a freezing of 
the advected spectra may be a good approximation, at least 



until 0.3 AU. If one assumes no further evolution, the double- 
power-law measured by Helios at 0.3 AU, may result from 
the combination of the parallel and perpendicular reduced 
spectra. At low frequency the parallel spectrum dominates 
with its characteristic slope -1, while at higher frequencies 
its steep spectral slope (-2) is masked by the more energetic 
perpendicular spectrum (with a -5/3 slope). Our theory 
might well be tested by the measurements of Solar Probe Plus. 
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